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THE CORRELATION BETWEEN MULTIPLICITIES OF 
CLOSED GEODESICS ON THE MODULAR SURFACE 

MANFRED PETER 



1. Introduction 



One of the connections between Number Theory and Mathematical Physics that 
I emerged in recent years is Arithmetical Quantum Chaos. On the physical side we have 
^ ' the problem of how the notion of chaos in classical dynamical systems can be transfered 
^ ' to quantum mechanical systems. Here Gutzwiller's trace formula is a useful quantitative 
. tool which connects the lengths of closed orbits in the classical picture with the energy 
eigenvalues in the quantum mechanical picture. Unfortunately it is an asymptotic rela- 
tion without rigorous error estimates. But in the special case of motion on a surface of 
constant negative curvature generated by a discontinuous group, Gutzwiller's trace for- 
! mula is Selberg's trace formula and is therefore exact. Furthermore, for particular choices 
^ I of groups connections with well understood number theoretical objects can be exploited 
CN ' which are not available for general groups. 

Numerical experiments led to surprising results: For a generic group the eigenvalue 
statistics of the Laplacian on the surface seem to be in accordance with the Gaussian 
O I Orthogonal Ensemble (see Bohigas, Giannoni, Schmit But for arithmetic groups it 
^ ■ is closer to Poisson distribution ([^, Aurich, Steiner see also 0). It seems that the 
d • high degeneracy of the length spectrum is responsible for the Poissonian behaviour of 
^ ! the eigenvalues. Note that a similar phenomenon was observed by Hejhal and Selberg 
j>I I for quaternion groups ([|], Theorems 17.1 and 18.8). They used the high degeneracy 
of the length spectrum to prove exceptionally large lower bounds for integral means of 
the remainder term in Weyl's asymptotic law for the eigenvalues. Luo and Sarnak 
', used the same phenomenon in their study of the "number variance" (the mean square for 
the remainder term in short intervals) for general arithmetic (not necessarily congruence) 
groups. They proved that for arithmetic groups, the length spectrum without multiplici- 
ties has at most linear growth and conjectured that this property characterizes arithmetic 
groups, a fact later proved by Schmutz |jT6| in the noncompact case. 



In an attempt to understand the Poissonian behaviour for the arithmetic group P : = 
SL(2, Z), Bogomolny, Leyvraz and Schmit calculated the two point correlation function 
for the eigenvalues. Their arguments are not mathematically rigorous and are given in 
two steps: First Selberg's trace formula along with heuristic arguments is used to reduce 
the pair correlation function of the eigenvalues to that of the lengths of closed orbits. 
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Second a heuristic version of the Hardy-Littlewood method is used to express the latter 
correlation function as an infinite product with easily calculable factors. 

Making the first step in |^ mathematically rigorous is extremely hard without a new 
methodological tool. Selberg's trace formula seems too weak for this purpose. In the 
present paper it will be shown how the second step can be made rigorous with an approach 
different from |Q. 

There is another way of looking at the result of this paper. Since the structure of the 
Selberg trace formula is the same as that of Weil's explicit formula in prime number theory 
(see [0]) one can look upon the lengths of closed geodesies as the logarithms of some sort 
of generalized primes. Thus the theorem below is an analogue of the Hardy-Littlewood 
twin prime conjecture. 

In order to state the main result, let EI be the complex upper half plane and F : = 
SL(2, Z). This group acts discontinuously on EI by Mobius transformations. Let F := 
be the Riemann surface generated by F. The closed geodesies on F are in one-to-one 
correspondence to conjugacy classes of primitive hyperbolic elements in F. For n G N, 
n > 2, let g{n) be the number of closed geodesies on F which correspond to primitive 
conjugacy classes with trace n. Set a{n) := g{n) logn/n. Then the function a has mean 
value 1. Set a(n) := a(n) — 1. 



Theorem. For r G Nq, the limit 



lim — a(n)a(n + r) 

N^oo N ^ \ I \ I 



2<n<N 



exists. Its value is given by 



p b>l 



where for a prime p > 2 and b > 2, we have 



Arip) 



(p2 - 1)^ 



1 



p 



E 



(n^ - 4)((n + r)2 - 4) 



P 



{p^ — ly p 



2 „36-4 



n mod p 

2/(1 r = 0(/) 

-2p*-\ r ^O(p^), r = 0(p*~i) 

(f)V(l-l/p), r = ±A{p^) 

-(— )V"\ r ^ ±4(/), r = ±4(/-i; 



> . 



0, 



otherwise 



Furthermore, for b > 6, we have 

1/1, r = 0(2) 



Ari2) 

A.(8) = 0, A'(16) 



9 1 -1, r = l(2) 
1 



9- 16 



1, r = 0(16) 
-1, r = 8(16) 
0, r ^ (8) 



1, r = (4) 
-1, r = 2(4) 
0, r = 1 (2) 



Ar{32) = 0, 



3 



AJ2^ 



9-2 



26-4 



0, 



r = (2^) 
r = 2^-1 (2^) 
r = ±(4 + 2^'^) 
r = ±(4 + 2''-2 
otherwise 



(2^) 
_26- 



(2'' 



This theorem was aheady stated in [Q and was made plausible on heuristic grounds. 
In the present paper, a different approach is used which exploits the connection of a{n) 
with a sum of class numbers of primitive binary quadratic forms (see (|2.1|) below). The 
main step is to show that a is almost periodic and to compute its Fourier coefficients. 
Then the theorem follows from Parseval's equation. To this end a method is used that has 

It might be of interest to note that almost periodic 
— already found apphcations to trace formulas for 



already been applied in |]T3| and |T^ 
functions — albeit on the real line 
integrable geodesic flows (see for an overview). 

The present paper is organized as follows: In section 2 the function a is reduced to 
class numbers. Section 3 contains a short review of almost periodic arithmetical functions. 
Furthermore, the function a is shown to be almost periodic by approximating it with a 
carefully choosen periodic function. This will considerably simplify the computation of 
the product representation of 7(r). In section 4 this Euler product is derived and its local 
factors are computed in section 5. 

Acknowledgement. I would like to express my sincere gratitude to Prof. Zeev Rud- 
nick for bringing this problem to my attention. I gained much from conversations with 
him and from the stimulating atmosphere which he and his co-organizers created at the 
DMV seminar "The Riemann Zeta Function and Random Matrix Theory" . 



2. Reduction to class numbers 

In order to state the connection with class numbers we need some notation. Proofs of 
the number theoretic facts in this section can be found for example in [0, |Q] and ]TT| . 

The letter d will always stand for a positive non-square discriminant (i.e. d = 0, 1 (4)). 
A primitive binary quadratic form is a polynomial / = ax^ + hxy + cy'^ with a, 6, c G Z and 
gcd(a, 6, c) = 1. Its discriminant is d = h'^ — Aac. Two such forms /j, i = 1,2, are called 

equivalent if there is a matrix ^ ^ G SL(2,Z) with fi{x,y) = f2{ax + l3y,'yx + 6y). A 

fundamental theorem in number theory now states that the number h{d) of equivalence 
classes of primitive binary quadratic forms with discriminant d is finite. This number is 
one of the important quantities in number theory since it appears in a surprising variety 
of situations. 

The Pellian equation -u^ — dv"^ = 4 has infinitely many solutions in integers {u,v). The 
solution {ud, Vd) with Ud, f d > and Ud minimal is called the fundamental solution since 
all the other solutions can be generated from it by a simple law. Set := {ud + Vd\fd)/2. 
The trace of is defined as ti{ed) := Ud- 

The one-to-one correspondence between primitive conjugacy classes in F and closed 
geodesies on F is as follows: Every primitive hyperbolic P G F has two real fixed points 
(one of them can be oo). The orthogonal circle in EI which ends in these fixed points 
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induces a closed geodesic on F of length I where cosh(Z/2) = | Tr P|. The possible values 
for / and their multiplicities are described in 

Proposition 2.1 (flSl, Corollary 1.5). The lengths of closed geodesies on F are the num- 
bers 21oged with multiplicities h{d). 

From this proposition and the preceding description it follows that 

logn 



ain) 



n 

d: tr{ej)=n 



E M^)- (2-1) 



Quantitative results about class numbers are often derived using Dirichlet's class number 
formula. First we must define Jacobi's character and Dirichlet L-series. Define Xd '■'^ ^ 
{0, ±1} to be completely multiplicative with 







= d {p) is 


solvable, p ){d 


Xdij)) ■= < 




x"^ = d ip) is 


insolvable 






p\d 








d=l{K) ] 




Xd{2) := < 




rf = 5(8) \ 


, Xd{-1): = 






c/ = 0(4) J 





1. 



Xd is a Dirichlet character modulo d (the quadratic reciprocity law is used to prove the 
d-periodicity). For 3fJs > 0, the series 



L{,s,Xd) ■■= XI 



Xd[.n) 



n" 

n>l 



is uniformly convergent and defines a holomorphic function. 

Proposition 2.2 (Dirichlet's class number formula). For a positive non-square discrim- 
inant d, we have 

h{d) \oged = VdL{l,Xd)- 

In the next section this formula will be used to prove the almost periodicity of a. In 
principle this could be done by writing 



^(1, Xd) = } ^ h error 



l<n<N 

Note that the sum on the right hand side is a periodic function of d and can be used to 
approximate L{l,Xd)- But this procedure has three severe drawbacks: First the approx- 
imation is not particularly good. Therefore we will use a smoothed version of the series 
representation of L(l, Xd) instead. Second we must approximate a sum of values L(l, Xd) 
with the condition tr(ed) = n. Breaking up this condition into easier summations, as must 
be done to process it further, would make things much more difficult. Third we want to 
compute the Fourier coefficients of a and show that they are multiplicative. This would 
be near to impossible with the above approach. 



5 



Instead an approximating periodic function is used which already incorporates some 
sort of multiphcativity. From the multiphcativity of Xd the Euler product 

pS 



follows. Thus it seems more reasonable to use a partial product of this representation 
than a partial sum of the series as approximating function. 

3. Almost periodicity 



The standard reference for almost periodic arithmetical functions is ^7j. Here the 
necessary material will be reviewed briefly. 

Let q > 1. For / : N — > C, define the seminorm 

1/9 



(hmsupl Yl l/HI') 'e[0,oo]. 

^ ^-"^ l<n<N ' 



f is called g-limit periodic if for every e > there is a periodic function h with ||/ — < e. 
The set of all g-limit periodic functions becomes a Banach space with norm || ■ \\g if 
functions /i,/2 with ||/i — f2\\q = are identified. If 1 < gi < g2 < oo, we have 
pi D D X)92 as sets (but they are endowed with different norms!). There is the more 
general notion of g-almost periodic function which will not be used in this paper (they are 
defined as above but with arbitrary trigonometric polynomials for h instead of periodic 
functions). For all / G "D^, the mean value 



^'(fy-^ti^ E /(«) 



l<n<N 



exists. The space is a Hilbert space with inner product 

{f,h):=M{fh), f,heV\ 

For u G M, define eu{n) := e^'^'™, G N. For all / G T)^, the Fourier coefficients 
f{u) := M{fe-u), u G M, exist. For m Q, we have f{u) = (this comes from the fact 
that / can be approximated by linear combinations of functions with u G Q, and that 
M(e„e_u) equals 1 if f — n G Z and otherwise; for almost periodic functions it is no 
longer true). In "D^, we have the canonical orthonormal base {ca/b}, where 1 < a < 6 and 
gcd(a, b) = 1. 

Limit periodic (and almost periodic) functions have a couple of nice properties. They 
can be added, multiplied and plugged into continuous functions and, under certain con- 
ditions, the result is again a limit periodic (almost periodic) function. They have mean 
values and limit distributions. Here we will use Parseval's equation. In Corollary |4.2| we 
will prove that a eV^ for all g > 1. As a side result in section 5, we get M{a) = a(0) = 1. 
Thus for r G Nq, we have a, a+j. G (where a+r{n) := a{n + r)), and Parseval's equation 
gives 



7(r) = M{aa+r) = (a, a+r) - 2M{a) + 1 



b>l l<a<b: {a,b)=l b>l l<a<b: {a,b)=l 

The last equation follows easily since a is real valued. If we can compute a and show that 
it is multiplicative our main theorem follows. 

The approximation of a by periodic functions is done in two steps. First we bring 
Dirichlet L-series into the picture. For n & N, n > 2, define 



V 

d,v>l: dv'^=rfi—'i 

where we must remember that d will always run through non-square positive discrimi- 
nants. 

Lemma 3.1. For q > 1, we have \\a — f3\\q = 0. 

Proof. For d fixed, the powers e^, / > 1, of the fundamental unit give all solutions {u,v) 
of the Pellian equation u"^ — dv'^ = 4 with m, f > 1 by way of the rule 

u + v\/d I 
2 

For every such solution we have v\fd = \/v? — A. Thus Proposition p.2| gives 

h{d) \oged 



d,l>l: tr{e')=n 



Since 



1 , ; 1 , n + \/n^ - 4 



it follows from (|2.1| ) that 



m = ^ ^ log !i±^^ y h{d) 

d,i>l:tr(ejj)=ri 

n logK- + v/^)^(^)^^/logn ^ 



\Jrfi - 4 log n V n 

(i>l,/>2:tr(ej;)=ra 

Since h{d) ^ \fd (this can be seen, e.g. from Proposition and the estimates loge^ > 
log|(l + v^) and L{\,Xd) ^ logc?; the latter follows easily by partial summation from 
the orthogonality relation for characters), it follows that for rf, Z > 1 with tr(e^) = n. 



h{d) « ^d « {e\r = f + < n^^. 

a{n) <^ logn 1 < logn r(n^ — 4) -C^ 



d,v>l: dv^=in?—A 
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for all e > 0; here r denotes the divisor function and r(m) was used. This implies 

\(3{n)-a{n)\ 



4 \-i/2 



A , log|(l + v/l-4/n2) ^ . logn 

IH -1 • H — > 

V logn / \ \ J\ ^ 



which proves the lemma. □ 

In the crucial second step (3 is approximated by a sum of partial products of the Euler 
product of L(l, Xd)- For n > 2, P > 2, define 

d,v>l: dv^=n^—4:,p\v=>p<P p<P 

Then /3(n) - (3p{n) = A^p\n) + A'p\n), where 

d,?;>l: dv^=n'^—i,p\v for some p>P 



2 /I „l ^».^D ^ P / 



and 

Ag^(n) := 

d,?j>l: dj;2=n2-4,p|t;=>p<P ^ p<P 

Fix g e N with g > 1 and choose q' > 1 with l/(2g) + l/g' = 1. The next lemma shows 
that A^p\n) is negligible in the 2g-norm as P — oo. 

Lemma 3.2. For P >2, we have 

i/g' 

\\2q 



^v>P 

Proof. Holder's inequality gives 

1/9' / \ 1/(29) 



|A?H|< 



^ d,v>l: ^ ^ > 1: dT;^ ^^2 _4 / 

p\v for some p>P p\v for some p>P 



For X > 1, this gives 

2<n<x ^v>P ^ 2<n<x,d,v>l:dv^=n'^-4 

The second sum on the right hand side is 

~ const ■ x 

d,i>l:tr(e^)<x 

as X — s> OO (see [0). This means that the values Xd) are constant in the mean when 
ordered according to the sizes of their fundamental units. Thus the lemma follows. □ 
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In order to estimate Ap (n) we must compare Xd) with a partial product of its Euler 
product. This is done by comparing both terms with a smoothed version of the Dirichlet 
series for L{l,Xd)- Let > 1. Then 



where 



Ei(i(i.x-j-5:^^-"") 

d,v ^ i>l ^ 



d,f i>l:p|i for some p>P 



Agi(«) 



d,-!; «>l:p|«=>p<P 



here the conditions on d and v are as in Ap''(n). The third term can be estimated easily. 
Lemma 3.3. For P > 2 and x,N > 1, we have 



2<n<x ^ Z>VJV;p|/^p<P 

Proof. Since |e~" — 1| ^ m for < m < 1, we see that for n > 2 the inner sum in Ap'^'*(n) 



IS 



E {k-""-i|« E 7+ E t4 

i>l:p|i=>p<P l>VN:p\l^p<P 1<1<VN 



< 5Z y + iV-i/2=:ci(P,iV). 

l>VN:p\l^p<P 

Holder's inequality now gives 

^ iAg,'(„)r'« ^ ( ^ i)\(p.iv) 

2<n<x 2<n<x ^ d,v>l: dv'2=n^-4 ^ 



2<n<x ^d,-u>l:dt)2=n2-4 ^ ^ d,-!;>l: dii2 = 

«Ci(P,iV)2'' ^ 1, 

2<n<x, d,v>l: n'^-dv'^=4 



where q' > 1 was used. The last sum is 

1 ~ const • X 

d,l> I: tr{e\)<x 

(see [O or [|15[), and the lemma follows. □ 



(2 2) 

So far we did not use the oscillation of the Jacobi character. For the estimation of Ap'jy 
we must take it into account. 



Lemma 3.4. For l,v ^ N and x > 3, we have 
where K{1) is the squarefree kernel of I and e > 0. 

Proof. See [IH, estimate (2.7). □ 



Lemma 3.5. For P > 2 and x, N > 1, we have 

2<n<x l>P^i ^ ' 

where nil) := Ei,,...,i,>i:i,...i,=i'^- 

Proof. Split the first sum in Ap'^^(n) into two sums depending on whether v < ^Jn or 
V > ^Jn. Thus Ap'^''(n) = Ap'^'^-'(n) + Ap'^'^'*(n). A trivial estimate gives 

^IT (-) « E - E 7 « N^rin^-A) 

— V — L \ ri 

d,v>l: dv^=n'^-4,v>^ l>l 

and thus 

J2 «, (logNf" J2 n^'^-"^^'^ «, (logiV)^'^ (3.2) 

2<n<x 2<n<x 

since g > 1. Holder's inequality gives 



d,v>l: dv'^=n'^—A, v<^/n 



I 

l>l:p\l for some p>P 



Xdjl) ^^i/N 



<i E E ( E 

ii,?;>l: ^1)2=712— 4 ci,?;>l: (i?;2=n2— 4, ti<y'n i>l:p|Z for some p>P 

Thus for X > 1, 

2q 



E |Ai?'^;»(«)r' « E E ( E ^ ^- 

2<n<x 2<n<x d,v'>l: dv'^=rfl—'i,v<y/x i>l: p|i for some p>P 



■l/N 



E E E xAh---i^)- 

ll,... ,l2q'-Pi\li ior some pi>P ^ l-CvK^/x 2<n<x,d>l: dv^=n^—A 



Applying Lemma p.4| to the innermost sum gives the estimate 

« E ] r^-M E j4«o ^ ^ 'fr^ ^ 

1>P'^1 1<V<^ ^ ' h,...,l2q>l 1<V<^ 
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which together with (p.2|) proves the lemma. □ 



In order to estimate Ap'^'*(n) we must show that the error 

/(ci,iV):=L(l,x.)-$:^e-'^ (3.3) 

i>i ^ 

which comes from smoothing the Dirichlet series expansion of L(l,Xrf), is small for large 
N . This is done by representing I{d, N) as an integral over a vertical line in the critical 
strip {0 < 3?s < 1} and using information about the location of the non-trivial zeros of 
L{s,Xd)- The Dirichlet series for L{s,Xd) is absolutely and uniformly convergent on the 
line = 2. Sterling's formula gives 

r(s) < I^SsI'^ e-"l^^l/2 
for ci < < C2, \Qs\ > 1, with some constant c = c(ci, C2). Using Mellin's formula 

— / T{s)y"'ds = e-\ 7/ > 0, 

we get 

-1. r(s - 1) L(., X.) AT-M. = ^ e-'/^. (3.4) 

On the other hand, for 1/2 < 3?s < 2, we have 

L{s,Xd)<^{d\s\Y/' (3.5) 

(this is easily seen by partial summation). Thus the line of integration in ( |3.4|) may be 
moved to the line = k with some 1/2 < /t < 1. Taking into account the pole of the 
integrand at s = 1, and using the residue theorem we get for (|3.4| ) the expression 

T{s- 1) L(s, Xd) N'"'ds + L(l, Xd) 

and thus 

1 pK+ioo 

lid, N) = -— / T{s - 1) L{s, Xd) N-'-'ds. (3.6) 

We must now considerably reduce the exponent 1/2 of d in (|3.5| ). In order to see the 
principle let us first assume the Generalized Riemann Hypothesis which says that for all 
Dirichlet characters x modulo q, the L-series L{s,x) has only zeros with real part 1/2 in 
the critical strip < 3?s < 1. From this the Generalized Lindelof Hypothesis follows and 
in particular for all d > 1, = k and e > 0, we have 

L{s,Xd) {d\s\Y. 

From (^) it follows that for > 1, 

I{d,N) d'N''-\ 
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This shows that 



d,v>l: dv^=n^—A 

and for P > 2 and a;, > 1, we get 



- J2 l^^Hn)]'" <^^,x'^'N"^^--'\ (3.7) 



X 

2<n<x 



Here it is important that the exponent of is negative. Taking to be a power of x with 
small exponent therefore lets the contribution of Ap'j^{n) vanish as a; ^ oo. The next 
lemma gives an estimate which for our purposes is as good as (p.7|) and can be proved 
unconditionally. 

Lemma 3.6. There are < k, /i < 1 such that for P > 2, x, N > 1 and e > 0, we have 
i E \A^^;i]in)\'' <^,x'N"^^-~'^+x''-'+'{\ogix'N)Y\ 

2<n<x 

Proof. Choose 1/2 < o"o < 1 with /i := 8(1 — cro)/o"o < 1- Choose < k, < 1. Define the 
rectangle 

R^:= {seC\ao<^s< 1, \Qs\ < hg^x}. 

If d < x^ and L{s, Xd) has no zeros in then a standard argument (see, for example, 
Titchmarsh |18[, Theorem 14.2) shows that for = k, \Qs\ < (loga;)^/2, we have 

log |L(s, Xd)\ < loglog(x + 2) (log(x + 2))('-")/('-"°) < elogx + c(e) 
with some constant c(e) > depending on e. Together with ( ^.51 ) and ( |3.(j| ) this gives 

J(rf,Ar)<, [ {\t\ + iy e-^'^'^/^x' N^-^dt 



+ [ Itl'^ e-"l*l/2 (^1^1)1/2 jyK-l^^ 



'|t|>(logx)2/2 

x'A^""\ (3.8) 

Next we must show that L{s, Xd) cannot have a zero in too often. From zero density 
estimates it follows that 

#{(?T,, V, d) \ 2 < n < X, d,v > I, n"^ — dv^ = 4, L{s, Xd) has a zero in R^^ -C x'^^" (3.9) 

(see [^, Lemma 4.11 or [|T1|, estimate (2.6)). A trivial estimation of ( p3.3[ ) gives 

I{d,N) <^log{dN). (3.10) 



3), (|3:i0|) and Hdlder's inequality give 



2<n<x 2<n<x ^ d,v>l: dv'^=n^ -4 ^ ^ d,v>l: dv'^=n'^~4 ^ 



2<ri<a;, d,v>l: 2<n<a;, d,v>l: dv'^—n^ —4 

L(s,X(i) has no zeros in Rx L{s,X(l) has a zero in Rx 
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which proves the lemma. □ 
Now the resuhs are collected. 
Proposition 3.7. For P > 2, we have 

Proof. For x > 1, choose N := x^^^^''\ Then Lemmas p.3|, |3.5| and p.6| show that 



2<n<x' ^ ,>a,i/(i69).p|;^p<p ^ ;>p29 V ; 

+ -nogx?'' + x-^/^+' + + x^-^+^(logx)2''. 

Since the series ^i>i.p|;=^p<p 1/^ converges, we have for P > 2 fixed 

lyp2q y > 



Together with Lemma this proves the proposition. □ 



4. An euler product 

In this section we exploit the particular construction of /3p by writing it as a product 
of functions each depending only on a single prime. 

For p a prime, 6 G No and n e Z, set Ipb{n) := 1 if = 4 (p^^) and, in case p = 2, if 
(n^ — 4)2"^'' is a discriminant (for p > 2 this is automatically fulfilled). Set Ipb{n) := 
otherwise. Define 

P(pM) ■=^-b{^ ~- X{n2-4)p-2b{p)^ Ipb{n), n>2. (4.1) 
Lemma 4.1. For P >2, we have jSp = Y[p<pP{p)- 

Proof. In the definition of l3p{n), write v = Y[p<pP^'' with bp G Nq. There is a discriminant 
d> with dv"^ = ra^ - 4 iff p'^^p\n'^ - 4 for all p < P and d := {n^ - 4)^"^ = 0, 1 (4). Since 
p2bp = 1 (4) for 2 < p < P, the last condition is equivalent to {n^ - 4)2''^^^ = 0, 1 (4). If 
these conditions are fulfilled, we have (n^ — 4)p~^^p = dr"^ with rp E N, p /r^, for p < P. 
Thus X(n2_4)p-26p (p) = Xd(yP) for p < P. This proves the lemma. □ 

Corollary 4.2. For q > 1, we have a G V^. In particular, limp^oo/5p = a with respect 
to the q-norm. 



13 



Proof. For g G N with g > 1 fixed and q' > 1 with l/(2g) + l/q' = 1 it follows from 



Lemma and Proposition |3.?| that 

||«-/3p||2,«C2(P)^/^' + C3(P)i/(2'^); 

here 



C2l 



v>P 



as P — s> oo since g' > 1. Furthermore, 



C3(P) := ^ - 



IK{1) 



as P oo, since 



T2,(/) 



alP' ■ a 



ECL \ — V 



k2e 



a>l b>l 



< oo. 



^IKil) 

1>1 a,b>l: a squarefree _ 

Thus limp^oo !!« — f^p\\2q = 0. For / : N ^ oo arbitrary and 1 < gi < g2 < oo, we have 
ll/llgi < Il/llq2 by Holder's inequality. Thus limp-^oo — (^p\\q = for all g > 1. 

Since the 6-th summand of /3(p) is p^^+^-periodic (22^'^'^-periodic in case p = 2) and 
the series representing /3(p) is uniformly convergent, the function /3(p) is uniformly limit 
periodic, i.e. G P"; here is the set of all functions which can be approximated to 
an arbitrary accuracy by periodic functions with respect to the supremum norm. Since 



"D" is closed under multiplication it follows from Lemma that (3p G for all P > 2. 
This gives a ^ V for all g > 1. □ 

Next the Fourier coefficients of a are computed in terms of the Fourier coefficients of the 
/3(p). In particular, this will show their multiplicativity. 

Lemma 4.3. (a) For all primes p, we have P{p){0) = 1. 

(b) Forb E N, a E Z, gcd(a, b) = 1, choose G Z for allp\b such that J2p\b^pP~"^'^^'' = 
ab^^ (1). Then 

"(^) =nAS(-^)- 

p\b ^ 

Proof. From Corollary |4.2| it follows that for arbitrary < e < 1 there is some P > 2 with 
\\a — /3p\\i < e and ordp6 = for all p > P. For all p < P it follows from ([4.1|) that there 
is some Ip > ordpb and coefficients Cp(a*) G C, 1 < a* < p'^, such that 



l<a*<p'p 



where || ■ \\u denotes the supremum norm and e' := eP ( maxp<p ||/?(p)||u + 1) • From 



Lemma 14. 11 it follows that 



/3p-n 



a;/ph 



P<P l<ai<ph 



< e. 
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Thus 



- E n 

l<a*<p'p (p<P) P<P 



< 2e. 



For / G we have |/| < Furthermore, Yl,p<pO*pP = ^(1) iff a* = 

^^pip-ordpb ^pip^ p < p. Therefore the orthogonahty relation for the exponential func- 
tion gives 

'a' 



a 



- n Cp(apP'''" 
p<p 



< 2e. 



Similarly, \P(p){apP 



This gives 



-ordpb) - Cp{app^p-°''^p^)\ <t'ioip<P and thus 



p<p 



n Cp(apP'''" 
p<p 



ord„ b\ 



< e. 



a 



p<p 



pordp b 



< 3e. 



In the next section we will compute and thereby show that P{p){0) = 1. This gives 
fa) and 



a 

"'6 



n'^(P)(pOrdp( 

p\b 



.p" 



< 3e. 



Since < e < 1 is arbitrary, (b) follows. 



□ 



From Corollary it follows that (|3.1| ) holds. Here the series on the right hand side is 
absolutely convergent (plug in r = 0). Thus Lemma [4.3| gives 



where for p prime and 6 G N, we define 



l<a<p^,p]/a 



Jliriar /p^ 



5. Computation of the local factors 

The last step is to calculate In particular, this will show that /3(p)(0) = 1 which is 
left over from the proof of Lemma [4.3|. For p prime, 6 G No, define 



/3(p,b) I 



n] 



1 \-i 
1 - - X(n2-4)p-26(p)j V(n), n>2. 



Then /5(p) = X]fe>oP~''Ap.'')' where the series is uniformly convergent. The calculation will 
only be done for p > 2. The case j9 = 2 is similar but somewhat more elaborate. Let 
6, c G No, a G Z, p /a. 
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Case 1: 2b < c — 1. Since /3(p,b) is p^^'^^-periodic, we have l3(^p^b){a/p'^) = 0. 
Case 2: b = c = 0. Then 

n mod p 

1 I 1,2 

The cardinahty of the first set is {p — 3)/2 and that of the second is {p — l)/2 (see for 

example the proof of Lemma 3.3 in [1^). Thus P(pfl){0) = 1 — 2/{p{p'^ — 1)). 
Case 3: 6 = 0, c = 1. Define 



S^{a) := Yl 



2nian/p 

nmodp:x„2_4{p)=±l 



Then 



n mod p 



Case 4: 26 > c - 1, 6 > 1. Then 

■''^ ■'^ ± nmodp2i'+i:n=±2(p2i') ^ 



Setting n = ±2 + mp gives 

it m mod p 

where (-) denotes the Legendre symbol. 
Case 4.1: 26 = c — 1. We have 

gT27rmm/p = 1 ^ gTSTriam/p ^ j _^ 1 j _ 1 
m modp: (m/p)=l mmodp ^ 



m\ 1 



1 \^ T27Tiam/pP 

2 ^ Vp; 2' 



Set ep := 1 if p = 1 (4) and := z otherwise. The last sum can be reduced to the 
Gaussian sum associated to the Legendre character which can be computed explicitely 
(see for example |P, Chapter 2). This gives for the above quantity the value 

Therefore 
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Case 4.2: 26 > c - 1. Then 

it m mod p 

U-l/p 2 1 + 2 y 

2 /47ra\p^+p + l 
cos ' ' 



p2b+l \ pC J p + 1 

Now we can calculate the Fourier coefficients 



6>0 

Case 1: c = 0. Then 

6>1 

Case 2: c = 1. Since 5'+(a) + S~{a) — — 2cos(47ra/p), we get 



^(0) = 1 - + y -l-g-cos(O) ^ = 1. 



2 /47ra\p^+p + l 
cos 



2p fA'Ka\ s-^ 1 ^1 , , 1 /n^-4\ 



P^ ~ - r , r \ ~ r ~ 

± n mod p 



Case 3: c > 2. Then 



0<6<(c-l)/2 



1 2 /47ra\ p- + /; + 1 

+ E 



cos 



&>(c-l)/2^ ^ I 



1 / ep((^)e^"*''/*'' + (^)e-4"^'^/p'), c odd 



(p2 - 1) p3c/2~2 I 2 cos (^) , c even 



Finally can be computed. 

Case 1: c—\. Then 



ni,n2modp / ^ ^ l<a<p— 1 



g27ri(ni-n2+r)/p 
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The innermost sum is p — 1 for ni — ^2 + r = (p) and —1 otherwise. Since 

/n^ — 4\ p — 3 p— 1 



n mod p 



the value for Ar{p) follows as given in the theorem. 
Case 2: c > 2. Then for p /a, we have 



^^"^ (pO I (p2 - l)p3c/2-2 h + ( p ) 



g-87rio/p" 



Thus 



(p2_l)2p3c-4 V Vp/ V p'^ 

l<a<p'',pjCa 

A short calculation gives the value in the theorem. □ 
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